Study of ground state phases for spin-1/2 Falicov-Kimball model on a triangular lattice 
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Abstract 

^ The spin-dependent Falicov-Kimball model (FKM) is studied on a triangular lattice using numerical diagonalization technique and 
Monte-Carlo simulation algorithm. Magnetic properties have been explored for different values of parameters; on-site Coulomb 
correlation U, exchange interaction J and filling of electrons. We have found that the ground state configurations exhibit long range 
^Neel order, ferromagnetism or a mixture of both as J is varied. The magnetic moments of itinerant {d) and localized (/) electrons 
are also studied. For the one-fourth filling case we found no magnetic moment from d- and /-electrons for U less than a critical 
value. 

T—I Keywords: A. Strongly correlated electron systems; C. Triangular lattice; D. Magnetic phase transitions; D. Exchange interactions 
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1. Introduction 

The correlated electron systems like cobaltates in Si, 
Gdl 2 li its doped variant Gd^Hx SSi, NaTi02 iE 
SE] , MgV 204 iflll] etc. have attracted great interest recently 
as they exhibit a number of remarkable cooperative phenom¬ 
ena such as valence and metal-insulator transition, charge, or¬ 
bital and spin/magnetic order, excitonic instability and possible 
non-fermi liquid states fl]. These are layered triangular lat¬ 
tice systems and are characterized by the presence of localized 
(denoted by /-) and itinerant (denoted by d-) electrons. The ge¬ 
ometrical frustration from the underlying triangular lattice cou¬ 
pled with strong quantum fluctuations give rise to a huge degen¬ 
eracy at low temperatures resulting in competing ground states 
close by in energy. Therefore, for these systems one would 
expect a fairly complex ground state magnetic phase diagram 
and the presence of soft local modes strongly coupled with the 
itinerant electrons. It has recently been proposed that these sys¬ 
tems may very well be described by different variants of the 
two-dimensional Falicov-Kimball model (FKM) Hi on the 
triangular lattice. 

Originally the FKM was proposed to describe the metal- 
insulator transition in mixed valence compounds lll2lll3n . Later 
the FKM was used to study the tendency of formation of charge 
density wave (CDW) order as well 1 13, 16, 17]. Recently 


many correlated systems ^ 25 j. In order to describe both 
the charge and magnetic orders in a unified way we use a gen¬ 
eralized FKM Hamiltonian lE^ that includes spin-dependent 
local interactions: 


we have studied the ground state and finite temperature prop¬ 
erties of the FKM and its different extensions on the triangular 
lattice liilii. We have reported several interesting results 
like various charge order, metal-insulator transitions and re¬ 
solved the issue of spontaneous symmetry breaking (SSB) iH 
in the ground state and explored the metal-insulator transi¬ 
tion at finite-temperature 21, 22] in the different regime 
of parameters. In all these studies the spin-degree of freedom 
was ignored and the interactions between electrons were spin- 
independent. Recent experimental results show that a charge 
order generally occurs with an attendant spin/magnetic order in 


H = 


+ r^ij)fcrdja- + {U J) flxfia-dl^dja- 

{ij)cr icr 

+ G ^ fj _^fi-a-d]^dia- + Uf ^ fia-fia-fl-o-fi-o' 

icr icr 

(1) 


here {ij) denotes the nearest neighbor {NN) lattice sites. The 
d- , dia- (T , ficr) are, respectively, the creation and annihilation 
operators for d- (/-) electrons with spin cr = {/ X} at the site 
i. First term is the band energy of the c/-electrons and fi is the 
chemical potential. The hopping parameter f(,y) = f for AA hop¬ 
ping and zero otherwise. The interaction between c/-electrons 
is neglected in FKM as usual. The second term is the on-site 
interaction between d and /-electrons of same spin with cou¬ 
pling strength {U - J) (where U is the usual spin-independent 
Coulomb term and J is the exchange interaction; the term fol¬ 
lows from Hund’s coupling). The third term is the on-site 
interaction U between d- and /-electrons of opposite spins. 
Here J basically represents the spin dependent local interac¬ 
tions between localized (/-) and itinerant (d-) electrons that 
stabilizes parallel over anti-parallel alignment between /- and 
c/-electrons. Inclusion of the exchange or Hund’s coupling term 
enables us to study the magnetic structure of the /-electrons 
and band magnetism of the ^f-electrons. Fourth term is on-site 
Coulomb repulsion U / between opposite /-spins while the last 
term is the spin-independent, dispersionless energy levelE f of 
the /-electrons. 

There are some theoretical results available for the spin- 
dependent FKM on a bipartite lattice 0261 l27l] . A few ground 
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state charge and magnetic configurations exist for certain fixed 
values of U and J. There is hardly any study available for the 
spin-dependent FKM on non-bipartite lattices. Therefore, in 
the present work we take up model systems that represent lay¬ 
ered materials with triangular lattice (hence geometrically frus¬ 
trated). Within second order perturbation theory, the spinless 
FKM with extended interactions can be shown to map to an 
effective Ising model with antiferrmagnetic interactions in the 
large U limit lll9n . The AFM coupling on triangular lattice is 
frustrated and leads to large degeneracies at low temperature. 
It turns out that this frustration is lifted |28, 2^ in the higher 
order perturbation in ^ J^. Therefore it would be quite inter¬ 
esting to see the role of spin degree of freedom of electrons on 
the ground state properties on such lattices with different values 
of parameters U and J. We study FKM at different range of in¬ 
teractions U and J for different electronic filling fractions on a 
triangular lattice. 


2. Methodology 

All the interactions in the Hamiltonian H (Eq.l) preserve lo¬ 
cal occupation and spin of the /-electrons, i.e. the (f-electrons 
traveling through the lattice change neither occupation numbers 
nor spins of the /-electrons. The local /-elctron occupation 
number hf/a- = is conserved as jia-, //j = 0 for all i and 

cr. This implies that ojicr - fi^fia- is a good quantum number 
taking values only 1 or 0 as the site i is occupied or unoccupied 
by an /-electron of spin cr, respectively. Following this local 
conservation, H can be rewritten as 

H — ^ "t" Uf ^ COio-OJi^-Q- 

(ij)cr icr 

+ Ef ^ COia- (2) 

icr 

where hij{{ua-]) = [ - tij + {(U - -H and 

{Wo-! is a chosen configuration of /-electrons of spin cr. 

The Hamiltonian H in Eq.2 shows that the /-electrons act 
as an external charge and spin-dependent potential or annealed 
disordered background for the non-interacting c/-electrons. This 
external potential of /-electrons can be “annealed” to find the 
minimum energy of the system. It is clear that there is inter-link 
between subsystems of /- and c/-electrons. This inter-link is re¬ 
sponsible for the long range ordered configurations and differ¬ 
ent charge and magnetic structures of /-electrons in the ground 
state. 

We set the scale of energy with t^ij) = 1. The value of 
is chosen such that the filling is (e.g. Nf + Nj = N 

is one-fourth case and Nf + Nj - 2N is half-filled case etc.), 
where Nf = (N^ + N/^), Nd - {Nd^ + Nd^) and N are the 
total number of /- electrons, d- electrons and sites respec¬ 
tively. Eor a lattice of N sites the H{{u)a-}) (given in Eq.2) 
is a 2N X 2N matrix for a fixed configuration {Wo-}. Eor one 
particular value of Nf{- -H A/^), we choose values of A/^ 
and Nf^ and their configuration {wf} = {wi|, ui 2 -r,..., wjv|} and 
{w^} = {wq, aj 2 i, ■■ ■, Wivi). Choosing the parameter U and J, 


the eigenvalues Aia-{i - 1... A) of /z({a;o-}) are calculated us¬ 
ing the numerical diagonalization technique on the triangular 
lattice of finite size N{- L?,L - 12) with periodic boundary 
conditions (PBC). 

The partition function of the system is written as, 

Z = ^ 7> (3^ 

i^a-] 


where the trace is taken over the (f-electrons, p - l/k^r. 
The trace is calculated from the eigenvalues T/o- of the matrix 
h{{(x>a-]) (first term in Eq.2). The partition function can, there¬ 
fore, be recast in the form. 



(4) 


Now, the thermodynamic quantities can be calculated as av¬ 
erages over various configurations {wo-l with statistical weight 
P({Wo-)) is given by 






where the corresponding free energy is given as, 

>[ 


FiM) = -- 
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(6) 


The ground state total internal energy ^({wo-)) is calculated as. 


Nj 

E{{Ua-}) = lim A({Wo.}) = V Aia-iitOa-]) 

T -^0 ^ 

i<T 

coiQ-oJi —fj- + Ef ^ tu/fT" (7) 

i(T io- 


Our aim is to find the unique ground state configuration (state 
with minimum total internal energy E({a>a.})) of /-electrons 
out of exponentially large possible configurations for a cho¬ 
sen Nf. In order to achieve this goal, we have used classical 
Monte Carlo simulation algorithm by annealing the static clas¬ 
sical variables {a>a-] ramping the temperature down from a high 
value to a very low value. Details of the method can be found 
in our earlier papers [18,19, 2^ 21, 22. 311. 


3. Results and discussion 

We have studied the variation of magnetic moment of 
c/-electrons (nid - ^ ^f-electrons are 

spread over all A number of sites) and magnetic moment of 
/-electrons {nif — ^ /-electrons are con¬ 

fined to only Nf number of sites) with number of (f-electrons 

Nf 

Nd at a fixed value of U, U / and J for nf - I (nf-j^ — 
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Figure 1: (Color online) Variation of magnetic moment of ^^-electrons mj and 
/-electrons m f with number of ^/-electrons for n/ = l,U = 5,Uf = \0 and 
for J = 5 and 3. 
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Figure 2: (Color online) (a) Up-spin and (b) down-spin r/-electron densities 
are shown on each site for 7 = 5, [/ = 5, f// = 10, «/ = 1 and Nj = 55. 
The color coding and radii of the circles indicate the 7-electron density profile. 
Triangle-up and triangle-down, filled by black and red colors coiTespond to the 
sites occupied by up-spin and down-spin /-electrons, respectively. 


—^). We have also studied the density of c/-electrons 
at each site for the above case. Fig.l shows the variation of 
magnetic moment of tf-electrons (mj) and /-electrons (;«/) with 
number of c/-electrons (Nd) for two different values of exchange 
correlation J i.e. J - 5 and 7 = 3 at a fixed value of on 
site coulomb repulsion t/ = 5 and Uf - 10. From Fig.l one 
can note that when Nd - 144 the ground state is Neel ordered 
anti-ferromagnetic (AFM) in nature. The reason for this can 
be understood in the following way. It is clear From Eq. 1 that 
for U - J, there is no repulsion between d- and /-electrons 
of same spins on the same site. The repulsion between d- and 
/-electrons of the opposite spins on the same site is U. There¬ 
fore it is energetically favorable that a c/-electron of same spin 
(as that of /-) occupies the site. For the FM arrangement of 
/-electrons, all li-electrons occupying the sites shall have FM 
arrangement themselves. Similarly for the AFM arrangement 
of /-electrons, 7-electrons occupying the sites shall have AFM 
arrangement. For FM arrangements of d- and /-electrons, there 


is no hopping possible for 7-electrons due to Pauli's exclu¬ 
sion principle, but for AFM arrangement of /- and 7-electrons, 
there will be finite hopping of 7-electrons between neighboring 
sites. This hopping reduces the kinetic energy of 7-electrons 
and hence total band energy of 7-electrons. Hence AFM ar¬ 
rangement of spins corresponds to minimum energy. Thus the 
ground state is AFM. In fact, it remains AFM up to Nd - 70 
for J - 5 and up to Nd - 65 for J - 3. The magnetic moment 
for /-electrons starts increasing for Nd <10 (for J - 5) and for 
Nd < 65 (for J - 3). Fully FM state is observed at Nd < 40 
for J - 5 and at Nd < 45 for 7 = 3, as now 7-electrons find a 
plenty of sites (where no 7-electrons are present) to hop to. By 
changing the number of 7-electrons, basically we are varying 
the doping and the extent of doping caused the phase transi¬ 
tion from Neel order AFM to FM via mixture of both state. 
These type of phase transition also observed using band struc¬ 
ture calculation for Gd^Hx, by doping hydrogen as reported by 
authors SSi. 

For the above mentioned parameter value t/ = 5 and J - 5 
and U = 5 and 7 = 3, we have seen /-electrons configura¬ 
tion and 7-electrons density at each site for different Nd- Fig.2 
shows the density of 7-electrons at each site for U - 5, J - 5 
for Nd - 55 (say). Table 1 and 2, respectively summaries the 
density of 7-electrons for two cases t/ = 5 and 7 = 5 (Fig.3) 
and U - 5 and 7 = 3 (not shown here) with FM and AFM 
arrangement of /-electrons. The density of 7-electrons at each 
site strongly depends upon the value of exchange correlation 7. 
Let us compare the density of 7-electrons at sites with up-spin 
and down-spin /-electrons for a fixed value of Nd (say Nd = 55) 
and for 7 = 5 and 7 = 3 (given in Table 1 and 2). We note from 
Table 1 that the density of 7-electrons at sites where 7- and /- 
electrons have same spins is large for 7 = 5 and less for 7 = 3. 
This is expected (as seen from Hamiltonian (Eq.l)), because 
for [7 = 7 = 5, there is no on-site repulsion between 7- and 
/-electrons of the same spins, but for [7 = 5 and 7 = 3 there is 
finite repulsion between 7- and /-electrons of the same spins. 
Also from Table 1 and 2, we note that density of 7-electrons 
with up-spin (n^^) at sites with EM arrangement of /-electrons 
is lesser than that at sites with AEM arrangement of /-electrons. 
This is so because in later case the 7-electrons may hop to ei¬ 
ther empty sites or to sites having down-spin /-electrons. This 
hopping reduces the kinetic energy of 7-electrons and hence to¬ 
tal band energy of 7-electrons. Eor the EM case, the 7-electrons 
can hop only to empty sites and the reduction in its total band 
energy is less. Same is true for density of 7-electrons with 
down-spin (n^^). 

We have also studied the ground state magnetic phases for 
up-spin and down-spin /-electrons, magnetic moments of 7- 
and /-electrons and the density of 7-electrons on each site for 
the range of values of parameters [7, U / and 7 for two cases (i) 
rif + rid - 1 (one-fourth filled case) and (ii) tif + rid - 2 (half 
filled case). We have chosen large value of Uf so that double 
occupancy of /-electrons is avoided. 

3.1. One-fourth filled case (nj + nd — 1).' 

In Eig.3 the ground state magnetic configurations of up-spin 
and down-spin /-electrons are shown for [7 = 5, [7/ = 10 


3 

























Table 1: The density of t/-electrons with FM aiTangement of /-electrons 


Density of 7-electrons 

Density of 7-electrons 

f/ = 5 

with up spin (n^^) 

with down spin (n^^) 

Nd = 55 

Sites with /y 

Sites with fi 

Sites with fy 

Sites with fi 


0.3 

0 

0.001-0.03 

0 

7 = 5 

0.27 

0 

0.058 - 0.072 

0 

7 = 3 


Table 2: The density of tf-electrons with AFM arrangement of /-electrons 


Density of 7-electrons 

Density of 7-electrons 

t/ = 5 

with up spin (n^/^) 

with down spin (n^^) 

10 

11 

Sites with /y 

Sites with fi 

Sites with fy 

Sites with fi 


0.37 

0.06 

0.06 

0.41 

7 = 5 

0.32 

0.1 

0.1 

0.3 

7 = 3 
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Figure 3: (Color online) The ground-state magnetic configurations of /- 
electrons for n/ = 5 , nrf = U = 5, Uf = 10 and for various values of 
J. Triangle-up and triangle-down, filled by black and red colors correspond to 
the sites occupied by up-spin and down-spin /-electrons, respectively. Open 
green circles con'espond to the unoccupied sites. 


and for different J values. The ground state configurations are 
observed to be long range ordered Neel ordered AFM arrange¬ 
ment of spins (Fig.3(a)) or complete FM arrangement (Fig.3(b)) 
or mixture of both arrangements (Fig.3(c). Here we note that 
U = 3.1 is the critical value of the on-site Coulomb correlation 
below which we do not get finite magnetic moment. Complete 
AFM phases are observed below U - 3.1. 

The variation of magnetic moment of c/-electrons {mj) and 
magnetic moment of /-electrons (nif) with exchange correla¬ 
tion J is shown in Fig.4(a) and Fig.4(Z7) for U - 3.1 and Uj -1 
and for U - 5 and Uf - 10 respectively. We note that in both 
cases the magnetic moment of d- and /-electrons increases with 
increasing J. Complete FM phase is observed at a particular 
value of J (e.g. J — 0.25 for U - 5 and Uf - 10) and it con¬ 
tinues up to some critical value of J (up to J — 2 for U — 5 and 
Uf - 10). With further increase in J mixed states are observed 
and finally AFM state is observed for larger values of J. 

Figs.5(a) and 5{b) show the density of up-spin ^f-electrons at 
a fixed value of U - 5, Uf - 10 and for J — 0 and J - 3.5 
respectively. When 7 = 0 the interaction between d- and /- 
electrons is same irrespective of their spins so the density of 
7-electrons at sites occupied by f-electrons are same, while it 
is maximum at unoccupied sites. With the increase in 7 value 
density of 7-electrons at sites where /-electrons of same spin 
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Figure 4: (Color online) Variation of magnetic moment of 7-electrons and 
/-electrons ntf with exchange correlation J for ny = j.nj = | at (a) f/ = 3.1, 
Uf = l and (b) U = 5,Uf = 10. 


are present increases and at empty sites it decreases, because as 
7 increases, the interaction (U - 7) between 7- and /-electrons 
of the same spins decreases. 

These results can be explained as we have three kinds of 
competing energies in the system namely the kinetic energy 
of 7-electrons and on-site interaction energies, ‘t/’ between 7- 
and /-electrons of opposite spins and ‘{U - 7)’ between 7- 
and /-electrons of same spins. For 7 = 0, the on-site inter¬ 
action energies between 7- and /-electrons of opposite spins 
and 7- and /-electrons of same spins will be the same. Hence 
the ground state configuration is AFM type as possible hop¬ 
ping of 7-electrons minimizes the energy of the system. It is 
clearly shown in the variation of 7-electron density at each site 
in Fig.5(a). As a result both the and nif are zero (as shown 
in Fig.4 and 7). For finite but small value of 7, the on-site inter¬ 
action energy between 7- and /-electrons of same spins will be 
smaller in comparison to the on-site interaction energy between 
7- and /-electrons of opposite spins. So few sites with FM 
arrangement of spin-up /-electrons will be occupied by down- 
spin 7-electrons and few by up-spin 7-electrons (as nid is very 
small). With this arrangement there is finite hopping possible 
for 7-electrons which reduces its kinetic energy and hence total 
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Figure 5: (Color online) Up-spin rf-electron densities are shown on each side for 
7 = 0 (a) and J = 3.5 (b), keeping other parameters same i.e. U = 5, Uf = 10, 
nf = ^ and nj = ^. The color coding and radii of the circles indicate the d- 
electron density profile. Triangle-up and triangle-down, filled by black and red 
colors correspond to the sites occupied by up-spin and down-spin /-electrons, 
respectively. Open green circles correspond to the unoccupied sites. 



Figure 7: (Color online) Variation of magnetic moment of if-electrons and 
/-electrons irif with exchange correlation J for (a) [/ = 2, Uf = 5 and (b) 
U = 5, Uf = 10 and for «/■ = 1 and rij = 1. 


energy of the system. Therefore the ground state configuration 
is either mixture of AFM and FM types or completely FM type. 
As a result both the and mf increase with increasing J. Fi¬ 
nally for large value of J (J ~ U), the ground state is AFM in 
nature. We have explained that in earlier section. 


3.2. Half-filled case (rif + 114 — 2): 


(a) J = 0.25, Nj = 96, N,. =48 (b) J = 1, = 132, = 12 (c) J = 2, = 72, N, = 72 
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Figure 6: (Color online) The ground-state magnetic configurations of /- 
electrons for U = 5, Uf = 10 and for various values of J with condition 
nf + Hd = 2. Triangle-up and triangle-down, filled by black and red colors 
correspond to the sites occupied by up-spin and down-spin /-electrons, respec¬ 
tively. 


Shown in Fig.6 are the ground state, magnetic configurations 
of up-spin and down-spin /-electrons for t/ = 5, t// = 10 and 
for different J values. At 7 = 0, again regular Neel ordered 
AFM structure is seen. With increasing J, mixed phase of FM 
and Neel type AFM is seen (Fig.6(a)). On futher increasing J 
the phase becomes fully FM and remains the same up to a value 
of J nearly equal to 0.75. Then the mixed phase comes back on 
increasing J at around 7=1 (Fig.6(b)) and finally the ground 
state becomes Neel ordered AFM again for 7 > 1.80 and con¬ 
tinues up to 7 = 5 (Fig.6(c)). Similar FM, AFM and mixed 
magnetic configurations of up-spin and down-spin /-electrons 
also observed for U — 2, Uj - 5 and for different 7 values (not 
shown here). 

Corresponding variation of magnetic moment of 7-electrons 
(mf) and of /-electrons (nif) with exchange correlation 7 is 
shown in Fig.7(a) and Fig.7(b) for U — 2 and U/ - 5 and for 
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Figure 8: Up-spin d-electron densities are shown on each side for 7 = 0 (a) 
and 7 = 0.5 (b), keeping other parameters same i.e. U = 5,Uf = 10, «/ = 1 and 
tid = 1. The color coding and radii of the circles indicate the (7-electron density 
profile. Triangle-up filled by black color corresponds to the sites occupied by 
up-spin /-electrons. 


U — 5 and U/ - 10 respectively. As we have already explained 
that the magnetic moments of d- and /-electrons increase with 
increasing 7 in both cases. Increase in magnetic moment of 
/-electrons is observed when one goes from Neel type AFM 
to complete FM phase through the intermediate mixed phase. 
Similarly, 1114 also increases with 7 but the increase starts after 
a certain value of 7. Here we observe that 1114 and nif increase 
initially with increasing 7 but after reaching a maximum the 
moments drop down at larger 7 values and finally both and 
nif become zero. We have found that for large value of U, mag¬ 
netization decreases sharply with 7 in comparison to that for 
small values of U. 

Figs.8(a) and 8(b) show the density of 7-electrons at a fixed 
value of U - 5, U/ - 10 and for 7 = 0 and 7 = 0.5 respec¬ 
tively. For 7 = 0 the density of 7-electrons at all sites is same 
irrespective of spins of /-electrons. As value of 7 increases, 
density of 7-electrons at sites where /-electrons of same spins 
are present increase as compared to at sites where /-electrons 
of opposite spins are present. This is expected because as 7 in¬ 
creases, the repulsion between 7- and /-electrons of the same 
spins decreases and hence they prefer to sit on the sites where 
/-electrons of the same spin are present. 

In conclusion, the ground state magnetic properties of two 


5 



































dimensional spin-1/2 FKM on a triangular lattice for different 
range of values of parameters like d- and /-electron fillings, 
on-site Coulomb correlation U and exchange correlation J are 
studied. We have found that the magnetic moments of d- and /- 
electrons depend strongly on the values of J and on the number 
of (f-electrons Nd- We have seen for one-fourth filling that there 
are no magnetic moments of d- and /-electrons for U less than 
3.1. At half-filling we have found that the magnetic moments of 
d- and /-electrons decrease sharply with J at larger U in com¬ 
parison to smaller values of U. At both fillings we note that the 
density of c/-electrons depends upon the value of exchange cor¬ 
relation J. Also, various charge and magnetic ordered phases of 
the localized /-electrons in the ground state have been observed 
at different values of J. The ground state configuration is ob¬ 
served to be long range ordered (either in some form of AFM 
arrangement of spins, complete FM arrangement or a mixture 
of both). The magnetic moments of d- and /-electrons start 
to increase sharply with increasing J and persist up to a larger 
value of J for one-fourth filling in comparison to one-half filling 
case (see Fig.4(b) and Fig.7(b)). There is rarely any calculation 
available for the spin-dependent FKM on a triangular lattice. 
Our results may motivate further studies of the magnetic prop¬ 
erties of frustrated systems of recent interest like Cobaltates, 
Gdl 2 , NaTi 02 and MgV 20 A,. 
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